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ABSTRACT 

We  propose  a  new  model  for  the  initiation  of  a  solar  coronal  mass  ejection  (CME).  The  model  agrees 
with  two  properties  of  CMEs  and  eruptive  flares  that  have  proved  to  be  very  difficult  to  explain  with 
previous  models:  (1)  very  low-lying  magnetic  field  lines,  down  to  the  photospheric  neutral  line,  can  open 
toward  infinity  during  an  eruption;  and  (2)  the  eruption  is  driven  solely  by  magnetic  free  energy  stored 
in  a  closed,  sheared  arcade.  Consequently,  the  magnetic  energy  of  the  closed  state  is  well  above  that  of 
the  posteruption  open  state.  The  key  new  feature  of  our  model  is  that  CMEs  occur  in  multipolar  topol¬ 
ogies  in  which  reconnection  between  a  sheared  arcade  and  neighboring  flux  systems  triggers  the  erup¬ 
tion.  In  this  “  magnetic  breakout  ”  model,  reconnection  removes  the  unsheared  field  above  the  low-lying, 
sheared  core  flux  near  the  neutral  line,  thereby  allowing  this  core  flux  to  burst  open.  We  present  numeri¬ 
cal  simulations  that  demonstrate  our  model  can  account  for  the  energy  requirements  for  CMEs.  We 
discuss  the  implication  of  the  model  for  CME/flare  prediction. 

Subject  headings:  Sun:  corona  —  Sun:  flares  —  Sun:  particle  emission 


1.  INTRODUCTION 

It  is  now  widely  recognized  that  coronal  mass  ejections 
(CMEs)  are  the  most  important  manifestation  of  solar 
activity  that  drives  the  space  weather  near  Earth  (Gosling 
1993,  1994).  CMEs  are  huge  ejected  plasmoids,  often  with 
masses  >1016  g  and  energies  >1032  ergs,  and  frequently 
subtending  more  than  60°  in  position  angle  (see,  e.g., 
Howard  et  al.  1985;  Hundhausen  1997).  Recently,  certain 
LASCO  coronagraph  observations  from  SOHO  have  been 
interpreted  as  evidence  that  CMEs  can  be  global  dis¬ 
ruptions  completely  circling  the  Sun  (Brueckner  1996; 
Howard  et  al.  1997).  In  addition  to  being  the  primary  cause 
of  major  geomagnetic  disturbances,  CMEs  are  also  a  funda¬ 
mental  mechanism  by  which  the  large-scale  corona  sheds 
flux  and,  hence,  may  play  a  central  role  in  the  solar  cycle. 
Therefore,  an  understanding  of  the  mechanism  for  CME 
initiation  has  long  been  a  primary  goal  of  solar  physics 
theory. 

There  are  a  number  of  basic  observational  constraints 
that  any  theoretical  model  for  CME  initiation  must  satisfy. 
Since  the  plasma  beta  in  the  nonflaring  corona  is  typically 
small,  P  ~  2 ,  gas  pressure  alone  cannot  be  the  driver. 
Early  models  for  CMEs  proposed  that  the  eruption  is 
driven  by  explosive  flare  heating,  but  it  is  now  known  that 
many  CMEs  occur  with  little  detectable  heating,  especially 
those  originating  from  high-latitude  quiet  regions.  It  has 
also  been  proposed  that  CMEs  may  be  due  to  magnetic 
buoyancy  effects  (see,  e.g.,  Low  1994;  Wu,  Gou,  &  Wang 
1995;  Wolfson  &  Dlamini  1997),  but  this  would  imply 
CMEs  should  be  associated  with  large  masses  of  falling 
material.  During  prominence  eruptions,  material  can  some¬ 
times  be  observed  to  fall  back  onto  the  chromosphere,  but 
CMEs  often  occur  with  negligible  prominence  activation 
and  very  little  evidence  for  downward  moving  plasma. 


Coronagraph  observations  usually  show  all  of  the  CME 
plasma  moving  outward,  in  which  case  buoyancy  is  unlikely 
to  be  the  driver.  These  considerations  have  led  most  investi¬ 
gators  to  conclude  that  the  energy  for  the  eruption  must  be 
stored  in  the  magnetic  field. 

The  general  topology  of  a  CME  appears  to  be  that  of  an 
erupting  arcade  overlying  a  photospheric  neutral  line  (see, 
e.g.,  Martin  &  McAllister  1997).  The  magnetic  field  near 
CME  neutral  lines  is  observed  to  run  almost  parallel  to  the 
line,  which  indicates  that  the  field  is  far  from  potential.  This 
magnetic  stress  is  commonly  referred  to  as  “  shear  ”  and  is 
most  likely  the  source  of  energy  for  the  eruption.  Note  that 
the  actual  type  of  arcade  can  range  from  low-lying  kilo- 
gauss  fields  associated  with  filament  eruptions  and  flares  in 
active  regions  to  weak  high-latitude  fields  associated  with 
polar  crown  filaments. 

Since  many  CMEs  are  accompanied  by  the  ejection  of 
prominences  or  filaments,  it  appears  that  the  innermost  flux 
of  the  erupting  arcade  opens  out  to  the  solar  wind.  The 
question  of  whether  all  the  arcade  flux  opens  or  whether 
some  inner  field  is  left  closed  is  a  critical  one  for  theoretical 
models  (Aly  1991;  Sturrock  1991;  Wolfson  &  Low  1992; 
Mikic  &  Linker  1994).  All  of  the  field  above  the  occulting 
disk  of  present  coronagraphs  (>1.1  R0)  is  observed  to 
open,  but  there  may  remain  a  considerable  amount  of 
closed  arcade  flux  below  this  height.  Since  low-lying  active- 
region  filaments  sitting  on  the  chromospheric  neutral  line 
are  observed  to  be  ejected,  we  believe  that  for  some  erup¬ 
tions  all  of  the  arcade  flux  down  to  the  chromospheric 
neutral  line  must  open.  Others  have  argued,  however,  that 
the  field  reconnects  back  down  to  a  closed  state  just  as  the 
filament  begins  to  lift,  so  that  the  innermost  flux  near  the 
chromosphere  never  actually  opens.  This  issue  remains  to 
be  resolved  observationally. 
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As  will  be  discussed  in  detail  below,  requiring  all  the 
arcade  flux  to  open  imposes  an  extreme  constraint  on  theo¬ 
ries  for  CMEs.  It  is  much  easier  to  explain  the  eruption  if  a 
significant  fraction  of  the  flux  remains  closed.  The  key 
feature  of  the  model  presented  in  this  paper  is  that  it  can 
produce  an  eruption  in  which  even  the  innermost  arcade 
flux  at  the  photospheric  neutral  line  can  open  to  infinity 
and,  yet,  the  eruption  is  purely  magnetically  driven. 

Another  observational  constraint  on  theoretical  models 
is  that  the  stressing  of  the  coronal  magnetic  field  is  slow 
compared  to  characteristic  timescales  in  the  corona.  Photo- 
spheric  driving  velocities  are  typically  of  order  1  km  s-1, 
whereas  characteristic  coronal  speeds  are  of  order  several 
hundred  km  s-1  for  sound  waves  and  several  thousand  km 
s-1  for  Alfvenic  motions.  Therefore,  to  a  good  approx¬ 
imation  the  energy  for  a  CME  is  pumped  into  the  corona 
quasi-statically. 

This  quasi-static  evolution  of  the  corona  is  controlled  by 
two  processes :  the  displacement  of  magnetic  footpoints  by 
photospheric  and  subphotospheric  flows,  and  the  emer¬ 
gence  and  submergence  of  magnetic  flux  through  the  photo¬ 
sphere.  It  is  not  clear  which  process  is  most  important  for 
producing  the  shear,  or  whether  both  play  significant  roles 
depending  on  the  type  of  solar  region.  Many  CMEs  orig¬ 
inate  in  coronal  helmet  streamers,  which  often  swell  for 
several  days  before  the  eruption  (Hundhausen  1997).  Such 
swelling  often  is  attributed  to  shearing  of  the  magnetic  foot- 
points  (see,  e.g.,  Klimchuk  1990),  although  the  emergence  of 
new  flux  also  could  produce  this  effect.  For  the  model  pro¬ 
posed  in  this  paper  it  is  irrelevant  whether  the  shear  is  due 
to  photospheric  flows  or  to  flux  emergence.  The  key 
requirement  is  that  the  stressing  is  slow.  In  the  simulations 
presented  below,  boundary  flows  were  used  to  generate  the 
shear,  but  we  expect  the  same  results  if  shear  emergence  is 
used.  The  major  advantage  to  our  boundary  conditions  is 
that  they  are  straightforward  to  understand  and  completely 
physical.  As  emphasized  by  a  number  of  authors  (see,  e.g., 
Aly  1984;  Klimchuk  &  Sturrock  1989;  Wolfson  &  Verma 
1991),  a  low-fi  system  like  the  corona  is  highly  sensitive  to 
boundary  conditions  at  the  base,  and  the  application  of 
inappropriate  conditions  can  easily  lead  to  erroneous  con¬ 
clusions  about  the  equilibrium  of  the  coronal  field.  By 
driving  the  shear  with  slow  footpoint  motions  we  are 
assured  that  any  eruption  that  occurs  is  physically  valid  and 
not  an  artifact  of  peculiar  boundary  conditions.  Note  that 
the  boundary  conditions  we  use  are  exactly  the  ones  used 
by  Aly  (1984),  Sturrock  (1989),  and  others  in  deriving  the 
energy  limit. 

The  final  constraint  on  CME  theories  is  that  the  corona 
has  no  upper  boundary.  Field  lines  are  allowed  to  expand 
to  infinity  if  it  is  energetically  favorable  for  them  to  do  so. 
This  seems  like  a  trivially  obvious  statement,  but  it  has 
profound  implications  for  the  possible  evolution  of  the 
magnetic  field. 

Taken  together,  the  constraints  of  an  open  corona  and  of 
quasi-static  stressing  make  it  very  difficult  to  understand 
how  violent  eruptions  can  occur  on  the  Sun.  Intuitively  one 
would  expect  the  field  to  respond  to  slow  photospheric 
driving  by  simply  expanding  outward  to  infinity  rather  than 
by  undergoing  a  violent  eruption.  This  intuition  is  con¬ 
firmed  by  a  number  of  numerical  simulations  that  have 
shown  that  sheared  2.5-dimensional  force-free  bipolar 
arcades  in  Cartesian  geometry  do  not  erupt  as  long  as 
physically  appropriate  boundary  conditions  are  assumed 


(Klimchuk  &  Sturrock  1989;  Finn  &  Chen  1990;  Wolfson 
&  Verma  1991).  There  is  also  no  evidence  for  eruption  in 
2.5-dimensional  bipolar  arcades  in  spherical  geometry 
(Roumeliotis,  Sturrock,  &  Antiochos  1994).  The  sheared 
field  simply  expands  outward,  which  pushes  any  overlying 
flux  outward  as  well  and  approaches  the  open  state  asymp¬ 
totically  with  increasing  shear  (Sturrock,  Roumeliotis,  & 
Antiochos  1995).  Eruptions  can  be  obtained  if  significant 
resistivity  is  present  (Biskamp  &  Welter  1989;  Steinolfson 
1991;  Inhester,  Birn,  &  Hesse  1992;  Mikic  &  Linker  1994), 
but  contrary  to  observations  the  field  opens  only  partially 
and  disconnected  plasmoids  form  in  the  corona. 

The  basic  problem  encountered  by  all  of  the  numerical 
models  to  date  is  that,  irrespective  of  the  magnitude  of  the 
shear,  the  closed  configurations  that  the  models  compute 
contain  less  magnetic  energy  than  the  fully  open  state.  In 
fact,  Aly  (1984)  and  Sturrock  (1991)  have  claimed,  on  the 
basis  of  mathematical  arguments,  that  the  maximum  energy 
state  of  any  force-free  field  is  the  fully  open  field  (Aly  1984, 
1991;  Sturrock  1991).  Although  this  claim  has  not  been 
completely  or  rigorously  proved,  it  is  in  agreement  with 
every  numerical  simulation  (including  the  ones  in  this 
paper),  and  hence  seems  very  likely  to  be  correct.  In  order  to 
obtain  a  CME,  however,  we  need  more  energy  in  the 
preemption  sheared  state  than  in  the  final  open  state,  since 
energy  is  needed  not  only  to  open  the  field  but  also  to 
accelerate  and  lift  the  plasma  against  gravity.  The  energy 
going  into  the  plasma  can  be  very  large,  of  order  1032  ergs 
or  more;  hence,  the  energy  in  the  preeruption  field  must  be 
well  above  the  energy  of  the  fully  open  state — a  result  which 
appears  to  be  in  direct  conflict  with  the  Aly-Sturrock  limit. 

2.  THE  THEORETICAL  MODEL 

We  propose  in  this  paper  a  CME  model  that  can  resolve 
the  apparent  contradiction  between  the  Aly-Sturrock 
energy  limit  and  CME  observations.  Let  us  first  clarify  why 
the  Aly-Sturrock  model,  in  fact,  requires  that  previous 
models  be  in  conflict  with  observations.  Since  a  number  of 
observations  indicate  that  CMEs  can  occur  over  extremely 
long  neutral  lines,  sometimes  apparently  circling  the  entire 
Sun  (Brueckner  1996),  and  since  most  CME  and/or  flare 
associated  activity  appears  to  occur  along  a  single  neutral 
line  on  the  photosphere,  models  for  CMEs  and  eruptive 
flares  generally  assume  a  2.5-dimensional  bipolar  magnetic 
geometry  with  either  translational  or  azimuthal  symmetry. 
Hence,  the  magnetic  system  consists  of  a  single  coronal 
arcade  (see,  e.g.,  Roumeliotis  et  al.  1994;  Mikic  &  Linker 
1994). 

The  important  point  is  that,  for  a  single  arcade  system, 
the  observational  requirement  that  the  innermost  flux  near 
the  neutral  line  open  up  requires  that  all  the  flux  in  the 
system  open;  but  this  is  forbidden  by  the  Aly-Sturrock  limit, 
at  least  for  a  purely  magnetically  driven  eruption,  because 
no  sheared  field  state  can  have  enough  energy  to  open  up  all 
the  flux.  The  only  way  to  obtain  an  eruption  in  these  models 
is  by  the  formation  of  closed  plasmoids  in  the  corona  so 
that  some  of  the  flux  is  disconnected  from  the  photosphere. 
The  Aly-Sturrock  arguments  do  not  apply  to  topologies 
with  disconnected  flux,  but  in  the  real  three-dimensional 
corona,  reconnection  will  not  result  in  the  formation  of 
disconnected  flux  from  the  photosphere.  In  three- 
dimensions,  the  Aly-Sturrock  energy  limit  will  apply  and  is 
likely  to  suppress  any  eruption  of  a  single  arcade  even  in  the 
presence  of  reconnection,  as  in  the  so-called  tether-cutting 
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model  (see,  e.g.,  Sturrock  1989;  Moore  &  Roumeliotis 
1992).  Furthermore,  it  is  very  difficult  to  see  how  three- 
dimensional  reconnection  in  a  single  arcade  can  lead  to  the 
opening  of  the  innermost  flux  near  the  neutral  line.  We 
conclude  that  any  single  arcade  model  is  doomed  to  failure. 

Let  us  consider  instead  a  multiflux  system  such  as  the 
field  shown  in  Figure  la.  In  this  case  we  have  three  neutral 
lines  on  the  photosphere  and  four  distinct  flux  systems:  a 
central  arcade  straddling  the  equator  ( blue  field  lines),  two 
arcades  associated  with  the  neutral  lines  at  ±45°  (green 
field  lines),  and  a  polar  flux  system  overlying  the  three 
arcades  ( red  lines).  Note  that  there  are  two  separatrix  sur¬ 
faces  that  define  the  boundaries  between  the  various  flux 
systems,  and  a  null  point  on  the  equatorial  plane  at  the 
intersection  of  the  separatrices. 

Suppose  that  as  a  result  of  large  shear  concentrated  at  the 
equatorial  neutral  line  a  CME  occurs  and  the  innermost 
flux  of  the  central  arcade  opens  up.  The  difference  now  is 
that  even  if  the  system  is  strictly  2.5-dimensional,  the  inner¬ 
most  flux  can  open  down  to  the  neutral  line  at  the  equator 


without  opening  all  the  flux  in  the  system  and  without  the 
formation  of  disconnected  flux.  Clearly  there  is  no  reason 
for  a  shear  that  occurs  only  near  the  equator  to  affect  the 
high-latitude  (green)  arcades;  therefore,  their  flux  should 
remain  closed  throughout  the  evolution.  This  is  not  signifi¬ 
cant  energetically,  however,  since  the  green  arcades  are  not 
expected  to  play  an  important  role  in  constraining  the  erup¬ 
tion.  The  key  point  is  that  much  of  the  flux  of  the  central 
(blue)  and  overlying  (red)  systems  also  can  remain  closed 
during  eruption.  If  reconnection  between  red  and  blue  field 
lines  occurs  at  the  null  point,  this  reconnection  will  result  in 
the  transfer  of  flux  from  the  blue  and  red  systems  to  the  side 
(green)  arcades.  Such  a  flux  transfer  allows  unsheared  blue 
and  red  flux  to  get  out  of  the  way  of  the  erupting  sheared 
core  flux,  but  to  still  remain  closed.  Note  that  this  process  is 
possible  only  in  a  multiflux  system.  In  a  single  arcade 
system  all  the  unsheared  flux  is  forced  to  open  because  it 
basically  has  nowhere  else  to  go. 

Flux  transfer  in  a  multipolar  system  also  allows  us  to 
understand  how  a  large  energy  excess  sufficient  to  power  a 


Fig.  1. — (a)  Initial  potential  magnetic  field.  The  field  is  symmetric  about  the  axis  of  rotation  and  the  equator,  so  only  one  quadrant  is  shown.  The 
photospheric  boundary  surface  is  indicated  by  the  light  gray  grid.  Magnetic  field  lines  are  colored  (red,  green,  or  blue)  according  to  their  flux  system.  Two 
types  of  blue  field  fines  are  indicated,  higher  lying  fight  blue  unsheared  field  and  low-lying  dark  blue  field  that  is  sheared  later  in  the  simulation.  ( b )  Force-free 
field  after  a  shear  of  71/8.  The  field  fines  shown  correspond  to  those  in  (a)  and  are  traced  from  the  same  footpoint  position  on  the  photosphere  as  in  (a),  (c)  As 
above,  but  for  a  shear  of  3it/8.  ( d )  As  above,  but  for  a  shear  of  tc/2. 
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CME  can  be  built  up  in  a  closed  sheared  field  and  yet  be 
compatible  with  the  Aly-Sturrock  limit.  The  key  point  is 
that  the  “fully”  open  state  is  not  unique  in  a  multipolar 
system.  To  illustrate,  let  £max  denote  the  energy  of  that  state 
in  which  all  the  flux  initially  in  the  central  and  overlying 
arcades  ( blue  and  red )  opens,  but  the  green  flux  remain 
closed.  When  the  shear  is  confined  only  to  the  arcade  at  the 
equator,  £max  is  the  appropriate  energy  limit  rather  than  the 
larger  energy  corresponding  to  all  the  flux,  including  green, 
being  open.  Note  that  the  state  with  energy  £max  can  be 
reached  by  a  purely  ideal  evolution  in  which  the  sheared 
core  field  expands  slowly  outward  pushing  all  the  overlying 
flux  out  to  infinity,  with  no  reconnection  at  the  null.  The 
side  flux  is  distorted  by  the  expanding  field,  but  it  does  not 
open.  The  final  state  has  a  current  sheet  at  the  equator 
reaching  down  to  the  photosphere  and  a  current  sheet  at 
the  top  of  each  side  arcade,  so  that  each  green  side  arcade 
becomes  a  helmet  streamer. 

We  expect  that  no  sheared  closed  configuration  can  have 
energy  above  £max,  so  if  it  were  the  only  state  accessible 
there  would  be  no  eruption,  just  as  in  the  single  arcade  case. 
If  we  allow  reconnection  at  the  null,  however,  then  other 
open  states  with  lower  energy  become  accessible.  For 
example,  if  the  shear  is  confined  to  a  small  latitudinal  band 
±  ©  near  the  equator  and  well  inside  the  separatrix  of  the 
blue  arcade,  then  consider  the  state  in  which  the  flux  inside 
±  ©  opens  while  as  much  as  possible  of  the  remaining  flux 
in  the  system  is  closed.  (It  may  not  be  possible  for  all  the 
unsheared  blue  flux  to  be  closed  because  there  may  not  be 
sufficient  red  flux  for  the  blue  to  reconnect  with.)  Let  this 
state,  which  has  the  minimal  amount  of  open  flux,  have 
energy  £min.  The  open  state  £min  looks  very  similar  to  the 
state  £max,  except  that  the  distribution  of  flux  between  the 
open  and  closed  systems  is  different.  We  can  consider  £max 
and  £min  to  be  the  maximally  and  minimally  open  states, 
respectively.  The  important  point  is  that  the  magnitude  of 
£min  can  be  much  less  than  £max,  because  the  state  £min  may 
contain  much  less  open  field. 

We  propose  that  the  energy  for  CMEs  and  eruptive  flares 
is  due  to  this  difference  between  £max  and  £min.  If  the  recon¬ 
nection  at  the  null  is  slow,  then  the  energy  in  the  sheared 
closed  configuration  can  rise  well  above  £min.  It  can  never 
go  above  £max,  but  depending  on  the  ratio  of  sheared  flux  to 
unsheared  flux  in  the  central  arcade  the  energy  of  the 
sheared  closed  field  could  be  much  larger  than  £min,  so  there 
could  be  a  large  amount  of  energy  available  for  eruption. 

In  this  model  the  unsheared  blue  and  red  fluxes  act  to 
confine  the  expansion  of  the  sheared  core  field.  We  expect 
that  a  CME  (or  eruptive  flare)  occurs  when  reconnection  at 
the  null  weakens  this  confinement  sufficiently  so  that  the 
sheared  field  starts  to  expand  outward  rapidly,  which  drives 
reconnection  ever  faster  at  the  null  and  “  breaking  out  ”  to 
infinity.  Note  that  this  magnetic  “  breakout  ”  model  natu¬ 
rally  implies  explosive-type  behavior. 

A  crucial  issue  for  our  model  is  the  rate  of  reconnection  at 
the  null.  If  this  rate  is  fast,  faster  than  the  rate  at  which  the 
system  is  being  sheared  (either  by  photospheric  flows  or  flux 
emergence),  then  the  magnetic  energy  always  stays  below 
£min  and  no  eruption  can  occur.  We  expect,  however,  that 
the  reconnection  will  be  slow  at  first  because  the  shear 
occurs  only  far  from  the  separatrices.  Under  these  circum¬ 
stances  a  quasi-static  evolution  should  be  possible  so  that 
only  weak  currents  develop  near  the  null  and  separatrices.  If 
so,  there  should  be  negligible  reconnection,  and  the  energy 


in  the  system  may  rise  well  above  £min.  In  the  next  section 
we  present  numerical  simulations  that  confirm  this  hypoth¬ 
esis. 

3.  NUMERICAL  SIMULATIONS 
3.1.  Force-free  Field  Calculations 
In  order  to  test  the  free  energy  issue  discussed  above,  we 
calculated  the  energy  of  a  sheared  multipolar  field  using  two 
procedures.  First,  we  calculated  the  energy  £min  by  deter¬ 
mining  the  minimum-energy  force-free  field  for  a  given  dis¬ 
tribution  of  shear  at  the  photosphere.  The  use  of  a  force-free 
calculation  needs  some  clarification.  Since  the  field  has  a 
null  point,  the  true  evolution  cannot  be  force-free  because 
any  finite  gas  pressure  will  dominate  near  the  null;  hence 
the  force-free  simulations  do  not  determine  a  physical  evo¬ 
lution  for  the  system.  We  use  the  force-free  code  only  as  a 
convenient  method  for  calculating  the  equilibrium  state 
corresponding  to  the  minimum  magnetic  energy  for  a  given 
shear  on  the  field  and  in  the  limit  of  vanishing  gas  pressure. 

The  important  feature  of  the  force-free  calculations  is  that 
they  include  the  effects  of  reconnection.  The  numerical  algo¬ 
rithm  that  the  code  uses  to  solve  the  force-free  equations  is 
basically  iterative  relaxation  (Yang,  Sturrock,  &  Antiochos 
1986).  Given  the  positions  of  all  the  field  line  footpoints  on 
the  boundary,  and  some  initial  guess  at  the  solution  in  the 
interior,  the  code  iterates  the  solution  monotonically 
toward  the  minimum  energy  state;  however,  the  code  uses  a 
multigrid  procedure  so  that  the  relaxation  is  first  performed 
on  very  coarse  grids  with  poor  resolution  and,  consequent¬ 
ly,  high  numerical  diffusion.  This  means  that  the  iteration 
process  allows  for  reconnection,  in  that  flux  will  transfer 
rapidly  from  one  flux  system  to  the  other  as  long  as  this 
transfer  is  consistent  with  the  boundary  conditions  and 
results  in  a  decrease  in  the  magnetic  energy.  Note  that  there 
is  no  real  diffusion  in  the  system,  however,  because  there  is 
no  slippage  of  the  footpoints  at  the  solar  surface.  The  shear 
is  strictly  enforced  by  the  boundary  conditions,  but  the  field 
is  free  to  reconnect  at  will  across  the  null  in  order  to  reach 
its  absolute  minimum  energy  state. 

The  code  solves  the  force-free  equations 

xB)xB  =  0  (1) 

using  an  Euler  potential  representation  for  the  field  (Yang  et 
al.  1986).  Azimuthal  symmetry  is  assumed,  so  that 

B  =  Va(r,  0)  x  V[0  -  y(r,  0)]  .  (2) 

The  advantage  of  the  Euler  representation  is  that  the  photo- 
spheric  shear  is  specified  directly  by  the  value  of  y  at  the 
photosphere  r  =  1.  By  fixing  y  there,  we  fix  the  shear  irre¬ 
spective  of  how  the  field  changes  in  the  interior. 

For  the  initial,  unstressed  (potential)  field  (Fig.  la),  we 
assumed  a  multipolar  flux  distribution  given  by 

m  sin2  9  (3  +  5  cos  29)  sin2  9 

a  r,  9)  = - + - — , - .  (3) 

r  2  r 

The  field  consists  of  a  dipolar  component  that  dominates  at 
large  r,  and  an  octopolar  component  at  the  surface.  In  the 
middle  corona,  r  ~  2,  the  field  appears  quadrupolar  with  an 
A- type  null  point  on  the  equatorial  plane.  Due  to  the  2.5- 
dimensional  symmetry  of  our  system,  the  normal  flux  at  the 
photosphere  is  unchanged  by  the  shear,  and  hence  the  value 
of  a  at  the  lower  boundary  is  constant:  a(l,  9)  =  (5/4) 
sin2  29.  Note  that  the  flux  maximum  a  =  1.25  occurs  at  the 
neutral  line  at  9  =  n/4. 
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From  equations  (2)  and  (3)  we  find  that 


Br(r,  9) 


1  8a 
r2  sin  9  89  ’ 


(4) 


so  that  the  normal  field  at  the  photosphere  is  given  by 
Br(l,  9)  =  10  cos  9  cos  29,  which  implies  that  there  are 
three  neutral  lines  located  at  9  =  n/4,  n/2,  and  3n/4.  The 
position  of  the  null  point  on  the  equatorial  plane  can  be 
determined  from  the  requirement  that  Bg  vanishes  there. 
From  equations  (2)  and  (3)  we  find  that  Be(r,  n/2)  —  r~3 
—  3 r~5,  and  consequently  the  null  is  located  at  r  =  31/2. 
The  value  of  the  flux  function  at  the  null  is  a(31/2,  n/2) 
=  2/[3  x  31/2]  =  0.385,  from  which  we  find  that  the 
separatrices  intersect  the  northern  hemisphere  at 
9  =  0.2941  and  1.277.  Note  that  we  have  chosen  the  form  of 
a  so  that  positive  flux  in  each  hemisphere  exactly  balances 
the  negative  flux  there,  while  a  vanishes  at  the  poles  and 
equator.  This  means  that  there  is  exactly  enough  red  flux  to 
reconnect  with  all  of  the  blue  flux,  if  such  reconnection  were 
energetically  favored. 

The  imposed  shear  is  assumed  to  be  antisymmetric  about 
the  equator  and  confined  to  a  small  latitudinal  band  there; 
i.e., 


n  m  =  /  xC(V  -  ©2)2  sin  \j/  ,  for  t/r  <  0  , 
n  ’  ’  \0  ,  for  ij/  >  0  ,  '  ' 

where  \j/  =  {n/2  —  9)  is  the  solar  latitude,  0  =  n/ 15  is  the 
assumed  latitudinal  extent  of  the  shear  layer,  and 
C  =  8.68252  x  103  is  a  normalization  constant  chosen  so 
that  y  =  x  at  the  latitude  of  maximum  shear,  \j/  =  0.094.  The 
form  of  y(l,  9)/x  is  shown  in  Figure  2.  Note  that  the  bound¬ 
ary  of  the  shear  layer  occurs  at  9  =  13n/30,  where 
a  =  0.207,  whereas  the  value  of  a  at  the  boundary  of  the 
blue  arcade  (the  separatrix)  is  0.385.  Hence,  only  about  half 
the  flux  of  the  central  arcade  is  sheared. 

We  solved  the  force-free  equations  above  in  the  domain 
(1  <  r  <  100,  0  <  9  <  n/2)  using  a  512  x  512  nonuniform 
grid.  Due  to  the  symmetry  of  the  system  about  the  equato¬ 
rial  plane,  only  one  hemisphere  needed  to  be  calculated. 
The  grid  was  selected  to  enhance  the  resolution  near  the 
solar  surface  and  near  the  equator;  in  particular,  the  grid 
points  were  spaced  uniformly  in  x  and  y,  where  x  =  In  r  and 
y  =  exp  (69/n).  The  boundary  conditions  at  the  pole,  0  =  0, 
and  the  equator,  9  =  n/2,  were  set  by  symmetry.  At  the 
inner  boundary  r  =  1,  the  shear  profile  y(l,  9)  was  specified 
using  the  form  above,  and  the  flux  distribution  a(l,  9)  was 
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Fig.  2. — Normalized  shear  profile  as  a  function  of  colatitude 


fixed  to  its  initial  value.  At  the  outer  boundary  r  =  100  both 
a(100,  9)  and  y(100,  9)  were  fixed  to  their  initial  value. 
This  meant  that  there  was  no  shear  at  the  outer  boundary, 
y(100,  9)  =  0,  which  is  a  good  assumption  because  the 
sheared  field  lines  were  initially  very  far  from  the  outer 
boundary.  However,  fixing  a  at  the  outer  boundary  means 
that  no  flux  is  allowed  to  exit  the  system.  This  is  not  valid 
since  as  the  inner  core  field  expands,  it  will  push  on  the 
outer  field  lines  so  that  even  the  flux  very  far  from  the  solar 
surface  will  move  outward.  The  amount  of  flux  at  r  =  100  is 
very  small,  however,  and  the  energy  associated  with  this  flux 
is  so  small  that  the  outer  boundary  conditions  have  a  negli¬ 
gible  effect  on  the  structure  and  energetics  of  the  system, 
even  for  shears  much  larger  than  the  values  we  discuss 
below. 

The  results  of  the  force-free  field  calculation  for  three 
shear  values,  x  =  n/S,  n/4,  and  n/2,  are  shown  in  Figures  lb, 
1  c,  and  Id.  The  field  lines  in  all  these  figures  were  traced 
beginning  at  exactly  the  same  footpoint  positions  as  in 
Figure  la,  in  which  there  are  four  footpoints  in  the  red 
system,  three  green  footpoints,  three  light  blue  (unsheared) 
footpoints,  and  three  dark  blue  (sheared)  footpoints.  The 
color  of  the  field  lines  in  Figures  1  b-ld  corresponds  to  the 
initial  color  of  the  footpoints,  irrespective  of  whether  recon¬ 
nection  has  occurred.  Even  for  a  low  shear  of  n/8  it  is 
evident  that  some  reconnection  has  taken  place  since  two  of 
the  red  and  two  of  the  light  blue  field  lines  have  now  joined 
the  green  system.  It  should  be  noted  that  although  the 
maximum  shear  is  n/8,  the  average  over  the  shear  layer  is 
only  about  half  this  value.  For  x  =  n/2,  almost  all  the  blue 
flux  has  reconnected  with  the  red  and  has  expanded 
outward  far  from  its  initial  position.  Therefore,  we  expect 
the  energy  of  the  n/2  sheared  state  to  be  close  to  that  of  the 
open  state,  Emin.  In  fact,  we  have  continued  the  force-free 
calculations  to  shears  twice  this  value  and  found  almost  no 
increase  in  the  magnetic  energy.  It  can  be  seen  in  Figure  3 
that  the  energy  of  the  force-free  field  appears  to  saturate  at  a 
maximum  value  of  approximately  6%  above  the  initial 
potential  field  energy;  i.e.,  Emin  =  1.06£pot.  Note  that  the 
6%  figure  refers  to  the  total  energy  of  the  initial  field,  which 
includes  all  the  unsheared  arcades.  If  we  consider  only  the 


Fig.  3. — Magnetic-free  energy  as  a  function  of  shear  for  the  force-free 
field  solution  ( lower  line)  and  the  MF1D  solution  (upper  line).  The  shear  is  in 
units  of  n  and  the  free  energy  is  expressed  as  percentage  of  the  initial 
potential  field  energy. 
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energy  initially  in  the  central  blue  arcade,  the  relative 
energy  increase  is  over  1  order  of  magnitude  larger. 

An  issue  that  requires  clarification  is  how  the  force-free 
calculations  determine  a  sequence  of  equilibrium  states 
leading  to  the  state  Emin.  The  code  finds  the  absolute 
minimum  energy  state  for  a  given  shear  at  the  base, 
assuming  that  reconnection  occurs  freely  at  the  null.  For 
small  shear  only  the  unsheared  light  blue  field  lines  recon¬ 
nect,  as  in  Figure  lb,  but  as  the  shear  increases  the  outer¬ 
most  sheared  field  line  eventually  reaches  the  null.  It  is 
important  to  emphasize  that  the  sheared  field  lines  ( dark 
blue  lines )  also  reconnect  freely.  This  is  evident  from  Figures 
lc  and  Id.  The  dark  blue  lines  have  reconnected  and  now 
have  one  of  their  footpoints  near  the  pole.  Although  these 
lines  have  reconnected,  they  have  not  lost  their  shear.  The 
footpoint  displacement  y  is  strictly  maintained  at  the  photo¬ 
sphere  by  the  boundary  conditions.  This  is  also  physically 
valid.  As  long  as  the  system  has  small  resistivity,  the  recon¬ 
nection  will  occur  in  a  region  of  negligible  volume  (a  current 
sheet),  and  the  field  reconnects  with  negligible  diffusion.  In 
this  case  the  reconnected  field  lines  maintain  their  footpoint 
displacement  (Karpen,  Antiochos,  &  DeVore  1996).  Note 
also  that  all  field  lines  with  a  footpoint  in  the  shear  zone, 
even  if  the  field  line  has  joined  the  green  system,  continue  to 
expand  outward  toward  the  open  configuration  as  the  shear 
increases.  In  the  limit  of  infinite  shear,  the  field  achieves  the 
open  state  Emin  in  which  all  of  the  field  lines  in  the  shear 
zone  (by  now  they  have  all  become  part  of  the  green  systems 
because  of  reconnection)  and  any  leftover  red  field  lines 
open. 

3.2.  MHD  Calculations 

The  results  of  our  force-free  calculations  determined  the 
magnitude  of  Emin.  The  key  question  is  whether  a  sheared 
coronal  field  will  evolve  to  a  state  with  energy  significantly 
above  this  value,  or  whether  reconnection  at  the  null  keeps 
the  energy  below  Emin.  In  order  to  calculate  the  physical 
evolution  for  the  field,  we  performed  a  fully  time-dependent 
simulation  using  a  2.5-dimensional  ideal  MHD  code  in 
spherical  coordinates  developed  by  DeVore  (1991).  The 
code  uses  a  multidimensional  flux-corrected  transport 
(FCT)  algorithm  for  transporting  the  conserved  quantities, 
and  it  maintains  the  divergence-free  condition  on  the  mag¬ 
netic  field  to  machine  accuracy. 

The  code  solves  the  standard  ideal  MHD  equations 
appropriate  for  the  solar  corona: 


|  +  v-  b»)  =  0. 

(6) 

|(P»)  +  v. 

'  (pvv)  +  VP  =  ^-(\  X  B)x  B-  pgQR% 

r 

?’ 

8U 

— +  V-  (Uv)  +  P\  ■  v  =  0  , 

ct 

(7) 

(8) 

and 

8B 

—  =  V  x  (v  x  B) , 

(9) 

where  p  is  mass  density,  v  is  velocity,  P  is  gas  pressure, 

B  is 

magnetic  induction,  gQ  is  the  solar  surface  gravity 
(2.75  x  104  cm  s-2),  Re  is  the  solar  radius  (7.0  x  1010  cm), 
and  U  is  the  internal  energy  ( U  =  3P/2). 


We  use  exactly  the  same  initial  magnetic  field  and  shear 
profile  as  in  the  force-free  calculations,  but  now  we  must 
also  specify  the  initial  plasma  distribution.  The  plasma  was 
taken  to  be  in  hydrostatic  equilibrium,  with  initial  tem¬ 
perature  and  density  given  by 

7  x  106  2  x  108 

T(r)  = - ■= —  K  and  n(r)  = - cm”3.  (10) 

r  r 

Since  the  field  drops  off  rapidly  with  radius,  these  forms  for 
T  and  n  were  chosen  in  order  to  keep  the  plasma  (1  = 
8nP/B2  from  becoming  too  large  at  large  distances  from  the 
surface.  Although  the  values  of  ft  in  our  simulation  are  not 
as  low  as  those  in  the  real  corona,  they  are  definitely  less 
than  unity  near  the  bottom  boundary.  The  values  for  T  and 
n  above  imply  an  initial  plasma  pressure  at  the  solar  surface 
of  5.5  x  10” 2  ergs  cm-3.  The  field  strength  at  the  surface 
ranges  from  a  value  of  B  =  Br  =  10  at  the  pole,  so  that 
/l  =  0.014  there,  to  B  =  Bg  =  —  2  at  the  equator,  where 
/l  =  0.35.  Once  we  begin  shearing,  the  value  of  /I  drops 
much  lower.  Of  course  the  system  is  always  high  beta  near 
the  null  point,  so  it  deviates  substantially  from  the  force-free 
case  there. 

We  use  the  same  nonuniform  grid  as  in  the  force-free 
calculations,  except  that  the  outer  boundary  is  placed  at 
r  =  10  rather  than  at  r  =  100.  As  will  be  shown  below,  the 
field  does  not  expand  outward  as  much  in  the  MHD  case, 
so  there  is  less  influence  from  the  outer  boundary.  Further¬ 
more,  we  assume  open  boundary  conditions  there,  so  that 
flux  and  mass  are  free  to  move  past  the  outer  boundary.  The 
use  of  a  smaller  domain  allows  us  to  have  higher  spatial 
resolution,  which  turns  out  to  be  the  limiting  factor  in  the 
MHD  simulations.  The  boundary  conditions  at  the  polar 
axis  and  the  equatorial  plane  are  determined,  as  before,  by 
symmetry  considerations.  At  the  inner  boundary,  the 
photosphere,  we  assume  a  line-tied  impenetrable  surface 
with  an  imposed  azimuthal  velocity  given  by  the  shear 
profile  above  and  with  a  sinusoidal  temporal  profile. 

A  critical  parameter  is  the  magnetic  Reynolds  number 
Rm.  Our  fully  two-dimensional  FCT  transport  algorithm 
uses  higher  order  differencing  so  that,  as  in  all  finite- 
difference  codes,  the  effective  Reynolds  number  is  quite 
high :  Rm  'p  10,000  as  long  as  the  magnetic  and  plasma 
structure  is  numerically  well  resolved.  Of  course,  these 
values  for  RM  are  still  small  compared  to  classical  solar 
values.  If  structure  develops  on  the  scale  of  the  grid  spacing, 
then  the  effective  RM  drops  to  become  of  order  the  number 
of  grid  points,  ~  500.  We  conclude,  therefore,  that  our  simu¬ 
lation  only  overestimates  (greatly)  the  true  rate  of  reconnec¬ 
tion  that  would  occur  on  the  Sun. 

We  drove  the  system  with  a  sequence  of  shear  motions 
applied  at  the  photospheric  boundary.  Each  shear  phase 
had  a  sinusoidal  time  profile  with  a  period  of  25,000  s  and  a 
maximum  footpoint  displacement  of  ti/8.  For  this  period 
and  amplitude  the  maximum  shear  velocity  at  the  photo¬ 
sphere  is  approximately  10  km  s_1,  whereas  the  Alfven 
speed  is  approximately  500  km  s_1  in  the  model  corona. 
Hence  the  evolution  was  nearly  quasi-static  until  the  very 
end  of  the  simulation  when  large  velocities  appeared  near 
the  null  because  of  reconnection. 

The  resulting  field  lines  at  the  end  of  the  first  shear  phase 
are  shown  in  Figure  4a.  These  field  lines  are  traced  from 
exactly  the  same  footpoint  positions  as  the  lines  in  Figure 
la.  We  note  that  the  field  lines  in  Figure  4a  have  pushed 
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Fig.  4. — MHD  solution  after  a  shear  of  (a)  ir/8,  ( b )  jt/4 ,  (c)  37i/8,  and  ( d )  n/2 .  The  field  lines  shown  are  the  same  as  those  in  Fig.  1. 


outward  as  a  result  of  the  shear,  but  unlike  the  correspond¬ 
ing  force-free  case  (Fig.  lb)  there  is  no  evidence  for  recon¬ 
nection  yet  in  the  MHD  results.  The  MHD  system  does 
appear  to  have  achieved  a  true  equilibrium.  We  let  the 
system  relax  for  an  additional  50,000  s,  and  the  resulting 
field  was  virtually  indistinguishable  from  that  in  Figure  4a. 
We  conclude  that  there  are  at  least  two  magnetostatic  equi¬ 
libria  that  the  system  can  achieve  for  a  given  shear.  One  is 
the  force-free  state  of  Figure  lb,  in  which  the  field  outside 
the  sheared  region  is  current-free  and  the  null  is  a  true 
A-point  with  right  angle  separatrix  lines.  This  state  has  the 
lowest  energy,  but  it  can  only  be  reached  with  reconnection. 
Another  solution  is  the  MHD  state  of  Figure  4a  in  which  no 
reconnection  occurs  at  the  null.  In  this  case  the  unsheared 
field  cannot  be  current-free  so  that  gas  pressure,  and 
perhaps  gravity,  must  play  a  role  in  the  force  balance.  We 
can  infer  from  Figure  4a  that  the  separatrix  surfaces  are  not 
quite  perpendicular  to  each  other,  which  indicates  there  is 
some  weak  current  near  the  null.  This  current,  however,  is 
less  than  1%  of  the  maximum  current  in  the  shear  region. 

Figures  4b,  4c,  and  4 d  show  the  results  of  the  MHD 
simulation  after  three  more  shear  phases.  The  solution  in 
Figure  4c  corresponds  to  a  total  maximum  shear  in  each 
hemisphere  of  3n/S.  Although  no  noticeable  reconnection 


has  occurred,  the  field  at  the  null  becomes  progressively 
distorted  from  a  right  angle  A- type  neutral  point.  The 
reason  for  this  is  straightforward.  To  relieve  their  stress, 
blue  field  lines  expand  outward  toward  an  open  configu¬ 
ration,  which  push  the  overlying  red  lines  outward  as  well. 
The  green  field  lines,  on  the  other  hand,  are  unstressed  and 
have  no  interest  in  being  dragged  out  to  infinity,  so  they 
simply  move  aside  by  pulling  away  from  the  neutral  point. 
Consequently  the  separatrix  lines  become  more  and  more 
oblique  as  blue  and  red  flux  push  toward  each  other  while 
green  pulls  away. 

We  can  gain  some  insight  into  the  evolution  of  the  struc¬ 
ture  near  the  null  point  by  considering  a  simple  analytic 
model.  Curvature  effects  can  be  neglected  sufficiently  near 
the  null,  and  a  locally  Cartesian  coordinate  system  can  be 
used.  Let  the  origin  of  this  coordinate  system  be  at  the  null, 
let  the  x-axis  be  horizontal  (parallel  to  the  equatorial  plane), 
and  let  the  y-axis  be  vertical.  If  the  field  near  the  null  has  no 
shear,  then  an  appropriate  form  for  the  potential  a  (see  eq. 
[2]) is 
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where  lx  and  ly  are  two  constants  that  determine  the  scale  of 
the  gradients  in  x  and  y,  respectively.  The  separatrix  lines 
are  given  by  a  =  0;  consequently,  if  lx  =  ly,  these  lines  are 
perpendicular  to  each  other  and,  as  we  show  below,  the 
current  vanishes.  From  this  form  for  a  we  derive  a  magnetic 
field, 

B"B°{t,*  +  T/)’  (12) 

and  an  electric  current  density, 

J  =  B0  ljyy*  z  ■  (13) 


force-free  and  MHD  solutions.  One  may  question  the  valid¬ 
ity  of  such  a  comparison  because  the  MHD  simulation  also 
includes  the  effects  of  the  plasma  energy;  however,  the 
change  in  plasma  thermal  energy  during  the  whole  shearing 
phase  (before  appreciable  reconnection)  was  less  than  10% 
of  the  magnetic  energy  change.  It  is  evident  from  Figure  3 
that  the  MHD  solution  energy  exceeds  Emin  even  for  a  shear 
of  7r/4,  while  the  free  energy  of  the  MHD  solution  is  approx¬ 
imately  twice  that  of  Emin  for  a  shear  of  3n/S.  Therefore, 
there  is  no  way  for  the  system  to  evolve  to  its  minimum 
energy  state  except  by  transferring  this  excess  energy  to  the 
plasma;  i.e.,  by  a  violent  eruption. 

4.  DISCUSSION 


Note  that  if  lx  =  /,„  the  current  vanishes  and  we  recover  the 
right  angle  X-type  null  point  of  a  potential  field.  If  the  scales 
are  not  equal,  there  is  a  finite  current  in  the  system  that 
implies  a  finite  Lorentz  force  given  by 


,14) 

Assuming  that  gravity  is  negligible,  the  Lorentz  force  must 
be  balanced  by  the  pressure  gradient,  which  implies  that  the 
pressure  is  given  by 


P  =  P, 


gp  ly  ~  K 
4; r  LL 


(15) 


where  P0  is  the  gas  pressure  at  the  null. 

Initially  the  field  is  current-free  and  lx  =  ly,  but  as  the  field 
is  distorted  by  the  outward  expanding  blue  flux  the  angle 
formed  by  the  separatrix  surfaces  becomes  more  and  more 
acute  so  that  lx  becomes  smaller  than  ly.  In  the  limit  lx  -4  ly, 
equations  (12)— (14)  above  imply  that  B  —>  (Bn  x/lx)y, 
J  -*•  (B0/lx)z,  and  P  -»  P0  —  (Bl/8n)(x2/ll).  In  other  words, 
the  null  deforms  to  a  classical  current  sheet.  This  trend  is 
clearly  evident  in  Figure  4. 

From  our  investigations  of  single  dipole  fields 
(Roumeliotis  et  al.  1994;  Sturrock  et  al.  1995),  we  expect 
that  at  large  shear  the  red  and  blue  field  lines  exponentially 
expand  outward  with  increasing  shear,  and  consequently 
the  current  sheet  thins  out  exponentially  fast.  This  implies 
that  reconnection  must  eventually  begin  to  occur  here,  and 
once  this  reconnection  starts  it  should  accelerate.  Our  simu¬ 
lation  shows  convincing  evidence  for  this  type  of  evolution. 
There  is  no  sign  of  reconnection  at  a  shear  of  3n/S  (Fig.  4c). 
We  let  the  field  relax  for  an  additional  30,000  s  after  this 
shearing  phase  and  saw  no  change  in  the  field  structure; 
hence  the  field  in  Figure  4c  is  also  a  true  equilibrium.  At  a 
shear  of  n/2,  however,  the  current  structure  at  the  null 
becomes  only  a  few  grid  cells  wide  and  reconnection  begins. 
For  example,  in  Figure  Ad  some  of  the  red  and  blue  lines 
clearly  have  reconnected  and  joined  the  green  flux  system. 
The  simulation  could  find  no  equilibrium  for  this  shear. 
When  we  tried  to  let  the  system  relax  by  continuing  the 
simulation  beyond  the  shearing  phase,  the  reconnection 
became  progressively  stronger  with  magnetic  islands 
appearing  at  the  interface  between  the  red  and  blue  systems. 
The  velocities  there  became  very  large  and  the  density 
plummeted,  halting  the  simulation. 

Although  our  code  could  not  simulate  the  actual  erup¬ 
tion,  it  did  prove  the  key  result  that  the  energy  of  the 
sheared  MHD  field  greatly  exceeds  Emin.  In  Figure  3  we 
compare  the  magnetic  energies  as  a  function  of  shear  for  the 


The  most  important  result  of  the  MHD  simulation  is 
that,  contrary  to  many  other  calculations,  current  sheets  do 
not  form  at  the  separatrix  until  late  in  the  shearing.  There 
are  two  reasons  for  this  result.  First,  no  shear  is  applied  at 
the  separatrix.  This  is  the  main  difference  between  this 
simulation  and  other  simulations,  including  our  own  (see, 
e.g.,  Karpen  et  al.  1996;  Karpen,  Antiochos,  &  DeVore 
1998).  If  a  shear  were  to  be  applied  at  the  separatrix  of  our 
configuration,  then  a  discontinuous  B,f>  would  be  created 
owing  to  the  discontinuity  in  relative  footpoint  positions, 
and  a  current  sheet  would  immediately  appear  just  as  in  the 
previous  simulations.  But  we  apply  a  shear  only  near  the 
neutral  line,  far  from  the  separatrix.  In  fact,  the  field  outside 
the  shear  region  should  always  be  nearly  current-free  for 
our  2.5-dimensional  system  since  a  quasi-static  force-free 
approximation  should  be  valid  everywhere  except  near  the 
null,  and  an  unsheared  2.5-dimensional  force-free  field  must 
be  potential.  This  is  exactly  what  our  simulation  finds ;  the 
current  is  completely  negligible  except  in  the  shear  region 
and  in  the  high-/!  region  around  the  null. 

Second,  discontinuities  do  not  appear  at  the  null  itself 
because  we  include  in  our  model  a  finite  plasma  pressure  so 
that  the  system  is  far  from  force-free  near  the  null.  In  order 
to  calculate  the  evolution  of  the  null  in  the  first  place,  the 
plasma  pressure  and  dynamics  must  be  fully  included. 
Without  the  plasma  the  null  could  not  even  change  its  posi¬ 
tion.  The  null  moves  outward  only  because  stress  is  trans¬ 
mitted  from  the  field  to  the  plasma  in  the  low-/!  region  and 
imparts  momentum  to  the  plasma,  which  then  moves  the 
field  in  the  high-/!  region  around  the  null.  The  plasma  pres¬ 
sure  is  also  essential  for  a  smooth  equilibrium  to  form.  As 
the  inner  flux  system  expands  outward  it  compresses  the 
plasma  in  the  high-/!  region  around  the  null,  which 
increases  the  plasma  pressure  there  until  force  balance  is 
achieved,  as  in  the  simple  analytic  model  described  above.  It 
is  only  when  the  null  region  is  squeezed  down  to  the  dissi¬ 
pation  length  scale,  which  corresponds  to  a  few  grid  points 
in  our  numerical  system,  that  significant  reconnection 
begins  and  equilibrium  becomes  impossible  to  maintain.  Of 
course,  on  the  Sun  the  magnetic  Reynolds  number  is  very 
large  and  the  dissipation  scale  is  many  orders  of  magnitude 
smaller  than  the  global  scale  of  the  magnetic  field. 

An  interesting  implication  of  this  discussion  is  that  the 
exact  energy  equation  for  the  plasma  must  play  a  crucial 
role  in  determining  when  eruption  will  occur.  In  the  simula¬ 
tion  we  assumed  a  simple  adiabatic  energy  equation  (8),  but 
radiation  losses,  thermal  conduction,  and  ohmic  heating 
may  all  be  important.  If  the  plasma  at  the  null  cools  as  it 
compresses,  the  width  of  the  null  region  can  decrease 
rapidly  to  the  dissipation  scale.  On  the  other  hand,  if  the 
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plasma  heats,  the  null  region  resists  compression  and  recon¬ 
nection  will  be  delayed.  Hence  we  have  the  very  interesting 
situation  that  the  dynamics  of  huge  phenomena  such  as 
CMEs  may  be  controlled  by  detailed  plasma  processes  that 
occur  in  relatively  tiny  regions. 

Another  interesting  property  of  our  “magnetic 
breakout  ”  model  for  CME  initiation/flare  triggering  is  that 
it  involves  several  observable  features  that  could  form  the 
basis  for  an  effective  prediction  scheme.  One  feature  is  that 
the  shear  is  concentrated  near  the  neutral  line,  which  is  in 
agreement  with  observations  (Schmieder  et  al.  1996)  and 
with  our  theory  for  prominence  formation  (Antiochos  et  al. 
1994).  Note  that  the  shear  responsible  for  eruptions  cannot 
be  due  to  differential  rotation  since  that  would  prefer¬ 
entially  shear  long,  high-lying  flux  rather  than  the  short, 
low-lying  flux  near  the  neutral  line.  Differential  rotation 
would  produce  the  opposite  shear  profile  from  that  shown 
in  Figure  2.  The  upcoming  Solar-B  mission  should  be  able 
to  measure  the  magnetic  shear  at  the  photosphere  with 
great  accuracy  and  to  determine  the  relation,  if  any, 
between  shear  distribution  and  eruptive  activity. 

Another  feature  of  the  model  is  that  multiple  flux  systems 
must  be  present  in  order  for  an  eruption  to  occur.  It  is  well 
known  that  large  eruptive  flares  only  occur  in  complex 
topologies;  (5 -spot  regions  in  particular.  It  can  be  shown 
that  a  (5-spot  has  a  three-dimensional  magnetic  topology  of 
a  four-flux  system,  exactly  analogous  to  the  2.5-dimensional 
model  discussed  here  (Antiochos  1998).  We  propose  that 
this  explains  why  (5-spots  are  so  flare  productive,  whereas 
bipolar  spot  regions  are  not.  Note  that  this  provides  a  clear 
distinction  between  our  breakout  model  and  the  usual 
tether-cutting  model  (see,  e.g.,  Sturrock  1989;  Moore  & 
Roumeliotis  1992).  Since  the  tether-cutting  model  involves 
only  a  single  arcade,  it  would  predict  that  a  strongly 
sheared  bipolar  spot  region  is  just  as  likely  to  flare  as  a 
(5-spot  region. 

CMEs  also  favor  magnetic  complexity.  Both  Mauna  Loa 
and  LASCO  observations  indicate  that  CMEs  are  more 
common  in  regions  where  the  coronal  field  has  a  multipolar 
structure  rather  than  a  simple  bipolar  structure  (McAllister, 
Hundhausen,  &  Burkepile  1995;  Schwenn  et  al.  1996).  In 
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addition,  observations  suggest  that  many  CMEs  involve 
more  than  one  neutral  line  and,  therefore,  more  than  one 
flux  region  (Webb  et  al.  1997). 

On  the  other  hand,  there  are  numerous  examples  of  high- 
latitude  CMEs  that  appear  to  involve  only  a  bipolar  mag¬ 
netic  topology.  It  would  seem  that  this  type  of  event 
contradicts  the  breakout  model;  however,  the  high-latitude 
CMEs  not  associated  with  active  regions  are  invariably 
slow,  with  height-time  profiles  similar  to  that  of  the  slow 
solar  wind  rather  than  flare-associated  ejection  (see,  e.g., 
Sheeley  et  al.  1997).  We  propose  that  this  type  of  CME 
should  not  be  considered  a  coronal  mass  ejection  but 
merely  a  coronal  mass  expansion.  The  key  point  is  that  if 
the  magnetic  field  expands  slowly  outward  because  of  shear 
while  the  plasma  maintains  its  temperature  due  to  coronal 
heating,  then  at  a  sufficiently  large  height  (a  few  solar  radii) 
the  plasma  will  begin  to  dominate  the  field  and  expand 
outward  indefinitely,  as  predicted  by  the  Parker  solar  wind 
model.  Most  CMEs  may  be  just  this  type  of  pressure-driven 
expansion,  but  obviously  this  process  cannot  explain  the 
fast  CMEs  that  appear  in  the  coronagraph  field  of  view  with 
velocities  much  faster  than  the  wind,  and  may  even  slow 
down  as  they  travel  upward  from  the  low  corona  (Sheeley  et 
al.  1997).  We  propose  that  fast  CMEs  must  be  due  to  mag¬ 
netic  breakout,  which  most  likely  occurs  in  strong  active 
region  fields. 

From  the  viewpoint  of  prediction,  the  most  important 
feature  of  our  model  is  the  reconnection  that  must  occur 
above  the  erupting  arcade.  This  reconnection  is  not 
expected  to  release  much  energy  since  most  of  the  free  mag¬ 
netic  energy  is  stored  in  the  low-lying  sheared  field.  Note 
also  that  the  reconnection  occurs  on  long  field  lines  far  from 
any  neutral  line.  Therefore,  it  is  unlikely  to  produce  signifi¬ 
cant  heating  or  strong  X-ray  emission,  but  it  may  be  detect¬ 
able  in  radio/microwave  emission  from  nonthermal 
particles  accelerated  by  this  reconnection.  We  believe  that 
radio  and  microwave  observations  will  provide  the  best  test 
of  the  magnetic  breakout  model  for  CME  initiation. 

This  work  has  been  supported  in  part  by  NASA  and 
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